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ABSTRACT 

Recently, Medvedev & Narayan (2001) discovered a new type of accretion flow, a hot settling flow 
around a rapidly rotating neutron star. The flow is cooling-dominated and energetically similar to the 
Shapiro, Lightman, & Eardley (1976, SLE) solution. Since the SLE solution is known to be thermally 
unstable, one might suspect that the new solution would also be unstable. However, due to the very 
high temperature of the accreting gas, thermal conduction is very strong and could suppress the thermal 
instability. We analyze the role of thermal conduction in both the hot settling flow and the SLE solution. 
In the hot settling flow collisions are very rare. Therefore, thermal transport occurs via free streaming 
of electrons along tangled magnetic field lines. We find that conduction is strong enough to make the 
flow marginally stable. In contrast, in the cooler SLE solution, conduction is via collisional, Spitzer-type 
transport. In this case, conduction is weaker, and we find that the SLE solution is thermally unstable 
even in the presence of conduction. 

Subject headings: accretion, accretion disks — instabilities — conduction — magnetic fields — stars: 
neutron — black hole physics 



1. INTRODUCTION 



Accretion flows around compact objects frequently ra 



The thermal stability of the Medvedev & Narayan 
(200l|) hot settling flow has not been addressed so far. 



diate significant levels of hard X-rays, indicating the pres- 
ence of hot optically-thin gas in these systems. A number 
of hot accretion solutions have been discussed in the lit- 



erature, e.g., the Shapiro, Lightman, & Eardley (1976) 



(SLE) sol ution, the advection-domina te d accretion flow 
(ADAF) ( [Narayan fc Yi 1994 )l995a| Jb|; [Abramowicrct 



al. 1995), and the convection-dominated accretion flow 
(CDAF) ( ^arayan et al. 2000| ; [Quataert fc Gruzinov 20001 ) 



These solutions are relevant for accretion onto a black hole. 



Recently, Medvedev fc Narayan (2001) discovered a so- 
lution that corresponds to accretion onto a neutron star 
(NS). We refer to this solution as a "hot settling flow," 
since the gas is hot (the temperature is nearly virial) and 
it "settles" onto the NS (the ratio of the radial velocity to 
the free-fall velocity decreases with decreasing radius). 

Not all hot accretion flows are stable. The cooling- 
dominatc d SLE solution has been shown to be thermally 
l unotab lo ( piran 1978 | ; | \^^andcl fc Liang 1991 | ; [Narayan fc Y\ - 



Energetically, this flow is very similar to the SLE solution, 
since the heat energy produced by viscous dissipation is 
radiated locally via bremsstrahlung. One might therefore 
expect the flow to be thermally unstable. However, this is 
not necessarily the case, as we show in this paper. 

The temperature of the gas in all these accretion solu- 
tions is very high, reaching almost the virial temperature 
in several cases. Thermal conduction is thus likely to be 
enormous and could have a signiflcant effect. The question 
that we try to answer in this paper is: How does conduc- 
tion affect the thermal stability of the hot settling flow and 
the SLE solution? 

In the hot settling flow, the mean-free-path is much 
larger than th e flow scale (the radius); therefore, the 
Spitzer theory (Spitzer 1962) of thermal conduction can- 
not be applied. Instead we need to apply the transport 
theory for coUisionless plasmas with tangled (but dynam- 



1005b [ ) and; honoo; unlilioly to oirict in naturo i Moro gon 
eralljr ^ it has been shown that an}" accretion flow in which 



ically unimportant) magnetic flelds (Rechcster fc Rosen- 
bluth 197SI; [Chandran fc Cowley 1998|; [Malyshkin fc Kul 



srud 2001[; [Malyshkin 200q). In tEecooler SLE solution 



heating balances cooling is thermally unstable if the cool- 
ing is due to brem sstrahlung emission (^hakura fc Sunyae'v 



1976[[Piranl978[ ). The ADAF solu tion, on the other hand 
is known to be thermally stable (Narayan fc Yi 1995b; 



Kato et al. 1996, 1997). In this solution, cooling is weak 
(ideally zero), and so the thermal energy of the flow is not 
radiated but is advected with the gas (hence the name). 
The CDAF is also believed to be stable, since in this flow 
again the thermal energy is advected by convective eddies 
and is either carried into the black hole or is radiated near 



the o uter boundary of the flow ( Ball, Narayan & Quataert 



2001) 



collisions are more frequent; hence it is appropriate to use 
collisional transport theory modifled for the presence of 
tangled flelds. 

The paper is organized as follows. We begin with a for- 
mal discussion of the thermal instability in §0, and discuss 
the stabilizing effect of thermal conduction in §|[ In SM 
and §H we study the effect of thermal conduction on the 
hot settling flow and the SLE flow, respectively, and we 
conclude with a discussion in §0. 

2. THERMAL INSTABILITY IN AN ACCRETION FLOW 



The physics of the thermal instability is simple (Field 



1965). Suppose a system is in thermal equilibrium, so that 
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the rates of heating and coohng per unit volume are equal: 
(5+ = Q~ . For simplicity let us take the heating and cool- 
ing rates to be functions of only the local temperature: 
Q+ oc T", Q- (xTP [a, (3 > Q for concreteness). 

Suppose, with increasing temperature, the cooling rate 
rises faster than the heating rate, i.e., (3 > a. Then a local 
perturbation which causes a small increase in the temper- 
ature will result in a net cooling of the gas: Q~ > Q^ . 
This will cause the temperature to return to its equilib- 
rium value, which means that the gas will be thermally 
stable. (It is easily seen that this is true also for a small 
decrease in the temperature.) On the other hand, if a > /3, 
the gas is thermally unstable. For instance, if the tem- 
perature decreases slightly, cooling becomes stronger than 
heating and the system deviates from its equilibrium in a 
run-away manner. 

To study the thermal stability of an accretion flow, 
we need to include additional physics, namely the ef- 
fects of shear and rotation. The shearing sheet ap- 



proxi mation (Goldreich fc Lyndcn-BcU 1965 ; Goldrcich & 
Treniidne 1978) is a convenient way of introducing the rele- 
vant physics without unnecessary technical complications. 
This approximation is quite accurate for perturbations on 
length scales much smaller than the local radius. 

Conventionally, the shearing sheet coordinates are 
Cartesian with x, y, z corresponding to the radial, az- 
imuthal, and vertical directions, respectively. These co- 
ordinates are appropriate for describing the motion of a 
parcel of gas whose geometrical size is small compared to 
the local radius, R, of the flow (i.e., x, y, z ^ R), so that 
the effects of geometry and curvature are insignificant. It 
is convenient to compare the wave- vector fc of a perturba- 
tion with 1/R and the frequency of a mode with the local 

Keplerian frequency i^K — \/GM/R^, where M is the 
mass of the central object. The shearing sheet approx- 
imation is accurate for "local" small-scale perturbations 
with kR ^ 1. Perturbations with kR ^ 1 are global; 
their properties may be understood only through a global 
stability analysis. 

We consider a shearing gas flow with unperturbed ve- 
locity given by 

Yo{x)^2Axy, (1) 

where 2 A ~ dVa/dx is the shear frequency and "hat" 
denotes a unit vector. Note that we have neglected the ra- 
dial velocity in the equilibrium flow since this component 
of the velocity is significantly smaller than the azimuthal 
velocity in both the hot settling flow and the SLE solution. 
To include the effect of rotation we assume that there is 
a Coriolis acceleration, described by an angular rotation 
frequency il = fi z. The vorticity and epicyclic frequency 
are then given by 



2B^2A + 2n, 



"I'Cp 



i = 2{nB) 



1/2 



(2) 



Both the hot settling flow and the SLE solution satisfy the 
Keplerian scaling, Q ex R^^^^. Therefore, for both solu- 
tions we have 2A = — (3/2)il, 2B — Vt/2 and Kepi = ^• 

We assume that perturbations in the flow have struc- 
ture only in the x direction, and we ignore motions in the 
z direction. We write the perturbations (represented by 
primes) in the velocity, density and sound speed as 

V'(x, t) — u{x, t) X + v{x, t) y, 



p{x,t) ^ poa{x,t), 
cl (x,t) = a^{x,t), 

where po ^^nd cj. are the equilibrium values of the density 
and the square of the sound speed. Note that we define 
Cs to be the isothermal sound speed, so that the pressure 
is written as p = pc^. By considering perturbations of 
the basic hydrodynamic equations, namely the continu- 
ity, radial momentum, azimuthal momentum and entropy 
equations, we obtain the following four equations. 



da 

dt 


du 
ox 


= 0, 
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+ cl 


dv 
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da da^ 
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PQCs 



,da 
'dt 



{Q+ + Q- 



(3a) 
(3b) 
(3c) 
(3d) 



where we have used d/dt — d/dt -f V^xd/dx ~ d/dt since 
the inflow velocity Vqx is set to zero in our approximation. 
Note that gravity does not enter the perturbed equations 
(it of course enters the unperturbed equation, where in 
equilibrium it cancels the centrifugal acceleration). In the 
shearing sheet approximation, the gravitational accelera- 
tion is assumed to be independent of x] therefore, it does 
not contribute to the perturbed equations. For simplicity, 
we have neglected viscosity in the azimuthal momentum 
equation, though we do include viscous dissipation in the 
energy equation through the terms Q^ and (5+ . 

For the heating and cooling rates, we make use of "re- 
alistic" expressions that represent the physics of viscous 
accretion flows. Thus we write 



Q+=a^ 



fdVf 



Qy 



PCs 



U,K V d^ 



AaA^^, Q- 



^ 



K 



-Cp' 



where a ~ 0.1 is the standard Shakura-Sunyaev vis- 
cosity parameter, V^y is the y-component of the unper- 
turbed velocity, and C is a constant. We leave the in- 
dex n in the cooling function unspecified for now, but we 
note that n = 1/2 corresponds to non-relativistic free- 
free (bremsstrahlung) cooling. In equilibrium, wc have 
Q^ + Qa = 0, and for the perturbations we find 



{Q+ + Q-y = -AA^a 



PoCs 



n 



K 



a 1 dv 
' c1 Adx 



We assume that the perturbations in equations (ph are 
of the form exjp[~iujt ~\- ikx). Substituting in the above 
equations and solving, we obtain the following dispersion 
relation: 



uj i{n — 1) 



7- 



'^cool 



■ UJ 



UJ ■ 



i{2B/A- 1)' 


k^cl 
k'cl 


'^cool 



where 



'^cool 






n 



K 



AA^a 9as2 



n 



K 



0, (5) 



(6) 



is the cooling (heating) time of the gas and s — VL/VIk is 
the dimensionless angular velocity of the gas. 

The dispersion relation (|5|) corresponds to purely ra- 
dial perturbations. The same relation can be used also 
for perturbations in the vertical direction, except that we 
must set Kepi = 0. Perturbations in the azimuthal di- 
rection are more complicated. Because of the shear, a 
non-axisymmetric wave packet is distorted as a function 
of time, and must be analyzed by special tec hniques which 
are beyond the scope of this paper (see, e.g., Toomre 1977; 
Goldreich fc Tremaine 19781 ). 

Equation (13) is a fourth-order polynomial and has four 
roots corresponding to four modes. A flow is unstable if 
any of the four modes grows with time, i.e. if the corre- 
sponding root has Im lo > 0. One of the roots of the dis- 
persion relation is always w = 0. This root corresponds to 
the viscous mode, which in the present case is particularly 
simple because we neglected viscosity in the momentum 
equation. It is easy to show that if we introduce viscos- 
ity into the momentum equation the viscous mode would 
become stable, i.e., we will obtain Im w < 0. We do not 
consider the viscous mode further in this paper. 

The physics of the remaining three modes may be under- 
stood by considering equation (0) in various limits. Con- 
sider first the limit fc — > 0. In this limit, two of the roots 
are given hy uj — ±Kcpi, corresponding to simple epicyclic 
oscillations. In the opposite limit fc — > oo, the same roots 
are given hy lu — ±7^'^Csfc, which shows that they cor- 
respond to sound waves. In the absence of heating and 
cooling (i.e. TcooI — > oo), we can obtain an exact solution 
for these roots which is valid for all fc: 



cpi 



7c; 



'.e 



(7) 



This is the standard dispersion relation for sound waves in 
a differentially rotating flow. The presence of 7 is because 
the relevant sound speed is the adiabatic sound speed, 
7^/^Cs (recall that Cg is defined to be the isothermal sound 
speed). 

The final root of the dispersion relation (0) corresponds 
to the thermal mode. In the limit fc ^ 0, we obtain 



LJ = Z (7 — 1) 



(1-n) 



We see that the mode is stable (for 7 > 1) 
unstable if n < 1. In the opposite limit fc -^. 

_ . (7 - 1) (2 - n - 2B/A) 



7 



'^cool 



(8a) 

if n > 1 and 
00, we find 

(8b) 



Now, the mode is stable if n > 2(1 — B/A)^ i.e. n > 8/3 
for our problem, and unstable otherwise. Note that an ac- 
cretion fiow that is cooled by free-free emission [n = 1/2) 
is unstable in both limits. 

The dashed lines in Fig. P show the real and imagi- 
nary parts of the various roots of the dispersion relation 
(||) as functions of fc for a realistic set of parameters. We 
have chosen K^pj = Q.b^\ which corresponds to a spin pa- 
rameter^ s = -\/0.5, and we have set the free-free cooling 
time to be Tcooi,ff — IO/JIk which corresponds to a ~ 0.1 

^The value of s 
higher than a typical spin of rotating neutron stars, which is around s 



(see equation B ) . Since we are interested in hot accretion 
flows with nearly virial temperature, we expect the sound 
speed to be comparable to the free-fall velocity; therefore, 
we have set Cg = ^IkR- Other parameters are 7 — 5/3 and 
n - 1/2. 

In Fig. m we do not show the trivial viscous root a; = 0. 
We label the two acoustic modes (propagating radially in 
opposite directions) as 1 and 2 and we label the thermal 
mode as 3. The frequencies of the acoustic modes are 
modified from the analytic solution given in equation (R), 
which was derived by neglecting cooling. The real parts 
of the roots are slightly perturbed, and both roots pick 
up an imaginary part. However, the imaginary parts are 
negative, which means that these modes are stable in the 
presence of cooling. The thermal mode, however, is un- 
stable (Im a; > 0) for all fc (Fig. ||), as expected from the 
asymptotic analysis presented earlier. We focus on this 
mode in the rest of the paper. 

3. STABILIZING EFFECT OF THERMAL CONDUCTION 

It is easy to see that thermal conduction will tend to 
reduce the thermal instability. An unstable thermal mode 
of wave- vector fc consists of a growing temperature pertur- 
bation of wave-length 27r/fc. Thermal conduction tends to 
smooth out this temperature perturbation through heat 
diffusion. If the rate at which the temperature perturba- 
tion grows is smaller than the rate at which it is smoothed 
out by conduction, then the instability will be suppressed 
and the mode will be stable. Otherwise, the mode will 
continue to grow, but at a somewhat reduced rate. 

The rate at which fluctuations are smoothed out by con- 
duction depends on the spatial scale of the perturbation. 
The smaller the scale (i.e. the larger the value of fc), the 
faster the conduction, and the greater the stabilizing effect. 
Thus, we expect conduction to stabilize thermal modes 
with fc greater than some critical fccrit- Our task in this 
section and succeeding sections is to estimate fccrit through 
a quantitative analysis. If we flnd that fccrit -R ^ 1, then 
we conclude that the flow is thermally unstable. On the 
other hand, if we find that fccrit^ < 1, we may reasonably 
claim that the fiow is thermally stable. Technically, for 
fc ~ 1/-R, we need to carry out a global analysis rather 
than the local analysis presented in this paper, but this is 
beyond the scope of the present paper. 

Let us write the heat flux q due to thermal conduction 
as 

gcond = -kVT, (9) 

where k is the thermal conductivity coefficient. Thermal 
conductivity in a dense, fully ionized gas is given by the 
Spitzer (1962| ) formula. 



Ksp « l.inkBVrX ^ 6.2 x IQ-^TJ^/^ erg/(s K cm). (10) 

Here vt = {kBTe/nie)'^/'^ is the electron thermal speed, T^ 
is the electron temperature (Tg = T for a one-temperature 
plasma), fcs is the Boltzmann constant, and 



A ~ IQ'^TJ/n cm 



(11) 



0.7 was chosen such as to illustrate the effect of rotation which is proportional to k^ . oc 



is the electron mean free path. Note that A is independent 
of the mass of the particle. 

This value of s is somewhat 
0.1. For such s the effect of rotation is negligible. 



In the coUisionless regime, i.e., when the mean free path 
of an electron becomes comparable to or larger than the 
temperature gradient scale A ^ Te/lVTej, equation ^ 
for the heat flux is no longer valid. For an unmagnetized 
plasma, the heat flux takes the following saturated form 
dCowie fc McKee 1977D, 



q,at ^ -Cpcl sgn(VT), 



(12) 



where C ~ 5 is a numerical constant whose exact value 
depends on the particle distribution function. This result 
is not relevant for our problem since our plasma is magne- 
tized. 

For a coUisionless magnetized plasma, thermal conduc- 
tion is anisotropic. Electrons stream freely along the fleld 
lines, and the parallel heat flux remains the same as for the 
unmagnetized case described above. However, the trans- 
verse heat flux is greatly reduced because electrons are 
tied to the field lines on the scale of the Larmor orbit. 
In fact, if the field is uniform and homogeneous, the per- 
pendicular thermal flux is identically equal to zero since 
electrons cannot move across the field lines. In a tangled 
field, however, electrons can jump from one field line to 
another and thus conduct heat perpendicular to the field. 
Since we are dealing with a turbulent accretion flow with 
a tangled magnetic fleld, this is the regime of interest to 
us. 

The ph ysics of this regime of conduct ion has been dis- 
cusse d bv Rechester fc Rosenbluth (1978[ ) and phandran fc 
Cowk y (1998 ), who identified two important efl'ects which 
are discussed in more detail in Appendix ^. 

First, since particles can move freely only along fleld 
lines, the characteristic effective mean free path is set by 
the correlation scale of the magnetic field Ib- In a hot 
accretion flow this scale is not known in general. How- 
ever, it is likely that turbulent motions in the fiow occur 
on a scale comparable to the local radius R, since this is 
the only characteristic scale in the problem. Very likely, 
the turbulent magnetic field will also have the same scale 
Ib ^ R- We parameterize this scale as Is = £,R- We ex- 
pect ^ < 1 because turbulent fluctuations cannot have a 
scale larger than the local radius of the flow. We assume 
<^ ^ 0.1 throughout the paper. 

Second, the magnetic fleld is inhomogeneous. Therefore, 
only a fraction z? < 1 of the particles will be able to pass 
though the magnetic mirrors that will be present in the 
fleld, and it is only these particles that transport energy 
beyond a distance '^ Ib- For magnetic fleld strength fluc- 
tuations 6B ~ (B) , the fraction of free streaming particles 
is estimated to be i? ^ 0.3. 

Typically, hot accretion flows are highly coUisionless, 
i.e., X ^ R ^ Ib- Therefore, we can write the thermal 
conduction coefficient as 



Kb — nkBVTlB ?? — 10 nkBVTR(,-i'&' 



(13) 



where ^_i = C/IQ-^ and d-i = -(9/10-^ Let us write the 
conductive heat fiux in a form similar to that used for the 



viscous stress, namely 



Qcoiid 



c^s del 
fix dx 



(14) 



where the dimensionless coefficient a^ is analogous to the 
Shakura-Sunyaev viscosity parameter a, and is given by 



R 
H 



^^F{e,p) ~ 10-2^_i^_i F(e,p). (15) 



Here we have used the fact that vt — Cse and H/R ^ 
Cs/vs ^ Cs/^kR, where H is the accretion disk scale 
height (in hot fiows, H ^ R) and vs is the free-fall speed. 
The quantity F{e,p) takes into account whether the con- 
duction is dominated by pr oton s or electrons; its numerical 
value is given in equation (A7). 

In the presence of conduction the energy equation has an 
additional contribution from the divergence of gcond . The 
perturbed energy equation (3d) then becomes modified to 



Po da^ 
j-l dt 



PoCs 



,d(T 
'di 



«3-^-«")'+'"i^-<''" 



This equation together with equations (3a-c) yields the 
following modified dispersion relation: 



2n2 



i{n-l) ik^R 



7-1 



"^cool 



'^cond 

i{2B/A - 1) 



Z Z ill 

uj - K^pi -k c. 



where 



Tcond = ^kR jOLc 



Pcl^O, (17) 



(18) 



is the conductive time scale. 

Figure |l| illustrates how thermal conduction modifies the 
thermal instability. The solid lines in the two panels show 
the real and imaginary parts of the roots of the disper- 
sion relation (|^ as functions of k. The parameters have 
the same values as before, namely Kepi/^A' = s = y/0.5, 
Cs/^kR = I7 Tcooi.e^K = 10, 7 = 5/3, n — 1/2. We have 
set the conductive time equal to rcooi,ff (for illustration), 
which corresponds to ac '^ 0.3. 

Figure ^ shows that the imaginary parts of all three 
modes decrease rapidly with increasing k- For the particu- 
lar parameters we have selected, the growth rate of the un- 
stable thermal mode (curve 3) goes to zero at /ccrit-R ^ 1-5, 
and the mode is stable for all k > /ccrit- 

We may also analyze equation (ll^) analytically. In the 
large-fc limit, the root corresponding to the thermal mode 
is equal to^ 



UJ 



(7-1) f{2^n-2B/A) k'^R^ 



1 



(19) 



'^cool '^cond , 

•^In deriving this equation from ( |l7| ) we have used the fact that the acoustic time-scale is, in general, shorter than the time-scale of the 
tiiermal mode, i.e., uj <C kcs, and wenave neglected K^pi as before. In this case we can neglect u)^ in the second brackets, so that equation 
(191) readily follows. It may seem that this procedure fails when the acoustic and thermal time-scales are comparable. This may happen when 
kR ~ 1 (for such perturbations, the sound crossing time is of order the dynamical time) and when r^ooi and T^ond a-r^ also comparable to the 



dynamical time, which is 
threshold itelf, Im oj = 0. 



o: 



Nevertheless, even in this case, equation lid) works fairly well near the stability threshold. Indeed, at the 



Since further the thermal mode frequpacy has no real part, we have oj ~ near the threshold; so we may safely 



neglect uj compared to k Cg in the second brackets of equation (hd) 



Clearly, for large k, conduction stabilizes the thermal 
mode, for the reasons explained at the beginning of this 
section. Using the above relation, we can estimate the 
critical fccrit above which all k are stable: 



•„2 i}2 'cond 



4 



13 Tc, 



6 Tcool ' 



(20) 



where we have substituted n ~ 1/2 (free-free cooling) and 
B/A = — 1/3 (Keplerian scaling). 

We should comment here that in the theory described 
above, the conductivity kb is a quantity averaged over 
many field correlation lengths. Therefore, the results of 
the stability analysis are valid only for perturbations on 
scales much larger than Ib- (This is somewhat inconsis- 
tent since we have assumed that Is = £,R ^ O.li?.) For 
small-scale perturbations with k 3> (^i?)~^, the local mag- 
netic field is nearly homogeneous. Therefore, thermal con- 
ductivity is anisotropic: its perpendicular component is of 
order of kb (equation [A3|), while the conductivity along 
the field lines is much larger: 



q^ _ Cpcl sgn(Vr) 
q± -ac{pc1/Q.K){dcl/dx) 



Cpcl 5 

IS. ^ ^1, 

acpcl{H/R) ac 



as follows from equations (h4h and (|12|). Thermal instabil- 
ity along the magnetic field lines is then suppressed more 
strongly than in the analysis presented above. 

4. THERMAL STABILITY OF THE HOT SETTLING FLOW 
SOLUTION 



The hot settling flow solution (Medvedev & Narayan 



2001) describes an optically thin two-temperature accre- 
tion flow onto a rotating neutron star. The flow extracts 
the rotational energy of the star and radiates it via free- 
free emission. The main properties of the solution are 
summarized in Appendix M. 

Because the hot settling flow is cooling-dominated and 
radiates by free-free emission, it is intrinsically thermally 
unstable, as follows from equation (pbh for 7 > 1, n = 1/2, 
and B/A = —1/3. However, thermal conductivity in the 
flow is enormous because of the very high temperature 
and large mean free path of particles. The coefficient of 
thermal conduction is estimated in Appendix M to be 



lO^^C- 



i^_iF(e,p), 



(21) 



where F{e,p) lies in the range 1 < F{e,p) < 15. Is this 
level of thermal conduction enough to stabilize the thermal 
mode? 



The relevant stability criterion is given in equation ( 2C 
However, before we apply this criterion, we need to allow 
for the fact that, as shown in Appendix pi thermal con- 
duction in the fiow is so strong that it modifies even the 
equilibrium structure of the flow. In particular, the cool- 
ing time (9) is modified as per the substitution given in 
equation (Bq) and becomes 



^cool 



(9as2 



T^^^ 



£"i'' 



where we have assumed that ac 2> as 
sonable for typical parameters, e.g. a ' 



(22) 
^ which is rea- 



0.1, 



ac ^ 0.1. F rom equation (18) and using the self-similar 



solution (Bl), we obtain the thermal conductive time'* 



3 

Oic 



'^cond — ^ ^}^ ; 



(23) 



Substituting TcooI and Tcond into the stability criterion I 
we then find 



fScrit^ 



26ac 



9(9as2 + 2ac) 




(24) 



that is, thermal modes with kR > 1 are stable. This sug- 
gests that the hot settling flow with thermal conduction is 
marginally stable to the thermal instability. 

Whether the mode kR = 1 itself is stable or not cannot 
be reliably determined from our local analysis. A global 
stability analysis is necessary to properly account for the 
effects of geometry and curvature, but this is beyond the 
scope of the present paper. 

5. THERMAL STABILITY OF THE SLE SOLUTION 



The SLE solution (Shapiro, Lightman, & Eardley 1976) 
was the first hot optically thin accretion solution discov- 
ered. This self-similar solution is best suited for accre- 
tion onto a black hole, though it could in principle be 
applied also for accretion onto a NS, by attaching a suit- 
able boundary layer at its inner edge. The SLE solution 
as originally envisaged by Shapiro et al. (1976) was cooled 
by Compton-scattering of soft photons from an o uter thin 
accretion dis k. A local form of the solution (e.g., Narayan 



& Yi 1995b) involves local cooling via free- free emission. 
Some properties of this solution are discus sed in Appe ndix 

SThe solution is thermally unstable (Piran 1978), as 
scussed in §2. The question we consider here is whether 
thermal conduction in the SLE solution is strong enough 
to eliminate the instability. 

In Appendix O we show that the SLE fiow is cooler than 
the hot settling flow and hence more coUisional. Therefore, 
the coUisionless thermal conductivity given by equation 
A3|) is not applicable. One should use instead equation 
To decide whether thermal conduction has any hope 
stabilizing the mode, we determine an upper boun d on 
the coefficient of thermal conduction. From equations (ICl) 
and (p7D we find 



KSLE 



<^ 1.1 X IQ^^dTit'/^C 



-5/2-15/8 



(25) 



The cooling time then follows from (p) for a Keplerian flow 

(5 = 1): 

(26) 
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Comparing equations (|9|) and (|lj), we write k = 
acpc^kB/i^K'iTT'p)- We then obtain a lower limit on the 
conductive time from equations (O) and (Ea): 
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* Alternatively, we recall that in the hot settling flow H/R • 



0-1, 

' 1 and H/R ~ Cs/vg ~ Cs/{\/2Q.kR). Then 



r^\ (27) 

"cond readily follows from equation 



Here we have used the valu es of plasma parameters for a 
two-temperature zone (C2), which yields a shorter Tcond 
than for a one-temperature zone. We notice that the con- 
duction time is significantly longer than the cooling time 
for typical parameters of the flow, e.g., to ^ 0.01, a ^ 
0.1, •& ~ 0.1, r <^ lO'*. Therefore, in contrast to the hot 
settling flow, thermal conduction in the SLE case is not 
strong enough to modify the equilibrium structure of the 
flow. We may thus directly apply the stability criterion 
([20|). We then find that the critical fccrit below which the 
thermal mode is unstable is given by 



kcritR ^ 180 (?9_im_2a-i) 



1/2^1/4 



(28) 



Since we find that kcritR is very large, we conclude that the 
SLE solution is thermally unstable even in the presence of 
thermal conduction. This is in contrast to the hot settling 
flow where conduction has quite a dramatic effect. The 
reason for the difference is that the SLE solution is cooler 
and more collisional and therefore has a significantly lower 
coefficient of thermal conduction. 

6. CONCLUSIONS 

In this paper we have investigated how thermal conduc- 
tion affects the stability properties of two hot, cooling- 
do rninMedj_^^tating__a£cretion_flows: the hot settling 
flow (iMedvedev fc Narayan 200l|) and the SLE slim disk 
(Shapiro, Lightman, & Eardley 1976). Without ther- 
mal transport, both flows are thermally unstable since 
they bal ance viscous heating with free- free brcmsstrahlung 
cooling (Piran 1978). Thermal conduction smears out 
temperature perturbations and can in principle suppress 
the instability. The stabilization effect is proportional to 
k'^B? , where k is the wave- vector of a perturbation and R 
is the radius, and is thus strongest for small-scale pertur- 
bations, kR 3> 1. 

To study the stability of the two rotating flows, we have 
used the shearing sheet approximation. We obtain the sta- 
bility criterion (EQ), which gives the critical wave- vector, 
fccrit J of a perturbation such that thermal modes with a 
larger spatial scale {k < fccrit) are unstable. We claim that 
an accretion flow is thermally stable if fccrit^ <, 1 and un- 
stable if kcritR ^ 1- Technically, for perturbations with 
k ~ \/R- the effects of geometry and curvature must be 



prope ply taken into account. This requires a global analy- 



sis, which we do not attempt. 

Because of the very high (nearly virial) temperature of 
the gas in the hot settling flow, the particle mean-free- 
path is much larger than the temperature gradient scale 
(which is comparable to R). In this regime, the Spitzer 
theory of thermal conduction cannot be applied. Instead, 
one should consider coUisionless transport theory includ- 
ing the effects of tangled (but dynamically unimportant) 
mag netic fields. The thermal fiux is then given by equation 
(|lj). We find that this flux is energetically very important 
in the hot settling flow and modifies even the equilibrium 
flow (see Appendix ^ . We allow for the modification and 
then analyse the thermal stability of the resulting fiow. 
We obtain the stability criterion ( p4| ) , which indicates that 
modes with kR ^ 1 are stable. This suggests that the hot 
settling fiow is thermally stable in the presence of thermal 
conduction. 

The SLE solution is cooler than the hot settling flow. 



and particle collisions are more frequent. The thermal 
flux in this flow is described by collisional transport the- 
ory suitably modified for the presence of tangled magnetic 
fields. We find that thermal transport in the SLE fiow is 
not very important and that the fiow is thermally unstable 
even when conduction is included (see eq ESh. 

We have thus answered the question with which we be- 
gan the paper (§1). However, in the process we have come 
up with an unexpected realization: conduction in hot di- 
lute accretion fiows may be so large as to substantially 
modify the equilibrium fiow. This raises three interesting 
issues. 

First, in the basic hot settling fiow without conduction, 
it has been shown that the energy source for viscous heat- 
ing of the gas is ultimately derived from the rotation of 
the central star (Medvedev & Narayan 2001). If, as we 



have shown here, conductive heating is more important 
than viscous heating, then one could ask where the energy 
for the conductive fiux comes from. Clearly, the energy 
must be derived from accretion. A complete global solu- 
tion, including conduction, should be able to trace the flow 
of energy from the hot inner regions of the accretion flow, 
where gravitational energy is released, to the hot outer 
regions, where the gas is approximately described by the 
self-similar settling flow solution. 

Second, if particle transport can have such a strong in- 
fluence on energy balance, what about angular momentum 
balance? Particles that move from one radius to another 
will certainly have some effect on the angular momentum 
of the gas, but the exact details are uncertain when parti- 
cle trajectories are constrained by the magnetic fleld. For 
simplicity, we have ignored in this paper angular momen- 
tum transport by the coUisionless particles. This effect 
should be included for consistency when the nature of the 
interactions is understood. 

Third, apart from the hot settling flow, two other hot 
accretion flows with nearly virial temperature are known: 
the ADAF and the CDAF. One suspects that energy con- 
duction (and angular momentum transport) by coUision- 
less particles might be important for these solutions as 
well, and might modify their equilibrium structures. We 
note that a CDAF is very different from an ADAF because 
convection has a strong effec t and changes the radial struc- 



ture of the flow dr astically ( Narayan et al. 200C ; Quataert 



& Gruzinov 2000). If we now add another form of trans- 



port, there could well be additional effects. 

This work was supported in part by NASA grant NAG5- 
10780. 



APPENDIX 

COLLISIONLESS THERMAL CONDUCTION IN TANGLED MAGNETIC FIELDS 

In this Appendix we discuss the essence of the thermal transport theory in coUisionless systems threaded by weak 
tangled magnetic fields. 

In the coUisionless regime, the mean-free-path A of an electron is comparable to or larger than the temperature gradient 
scale A > Te/lVTej. The heat flux in this case is ~ pc^. In a coUisionless magnetized plasma, electrons stream freely 
along the field lines only. In the transverse direction, they are tied to the field lines on the Larmor scale. In a tangled 
field, however, electrons can jump from one field line to another and thus conduct he at across the field. This is the re gime 



of interest to us. The phy sics of this regime of conduction has been discussed by Rechester fc Rosenbluth (1978 ) and 



Chandran fc Cowley (1998) 



Let us consider a tangled magnetic field with a characteristic correlation scale Is- We assume that the electron Larmor 
radius pe in this field is much smaller than the mean free path A. If / is the path length measured along a field line, the 
separation between two closely neighboring field lines grows as 

d(0~(i(0)e'/^^-, (Al) 

where d(0) is their initial separation and Lk is the Kolmogorov-Lyapunov length which, in general, depends on the field 
spectrum. Since there is only one characteristic scale in the problem, Ib, we expect Lk ^ Ib- Let us assume that at an 
initial instant, an electron drifted a distance d(0) ~ pe from its initial field line in the perpendicular dir ecti on. Once it 



starts moving along a new field line, its new path diverges from its initial path according to equation (Al). After the 
particle travels the distance 

Lrr^IbHIb/pc) (A2) 

measured along the field line, the separation becomes of order Ib- Any subsequent motion of the electron then becomes 
uncorrelated with its initial field line. The distance Lfm is called the Rechester-Rosenbluth length. 

If vt is the characteristic electron thermal speed, then from dimensional considerations it follows that the thermal 
flux is describ ed by equation IM) with an effective conductivity kb ~ nkBVTls- A more rigorous analysis ( Chandran &| 
Cowley 1998) shows that for large-scale perturbations, Mb < 1, one can take the ensemble average. Then the thermal 



conductivity becomes isotropic (assuming isotropic magnetic turbulence) and reads as 

min(A,ii?fl) 
Kb — nkBVTiB'ff 7 , (A3) 

where the function min(a, b) denotes the smaller of a and b. The quantity •& < 1 takes into account the fact that only a 
fraction of particles (those which are not trapped by magnetic mirrors) can move far enough to transport heat; the rest 
will remain trapped in their local magnetic wells and be unable to transport energy beyond a distance '^ Ib ■ We calculate 
^ below. 

There are three substantially different regimes of thermal conduction. First, in the coUisionless limit, A ^ Ljir, we have 
Kb ~ nkBVTlB'&i i-e., Ib plays the role of an effective mean-free-path and only non-trapped particles contribute. Second, 
in the semi-collisionless limit, pe ^ A < Lrr, we have kb ~ kspI^Ib/Lrr. That is, the standard Spitzer coUisional 
transport is affected by the presence of tangled flelds: the mean- free-path is set by collisions but most of the particles are 
still trapped in magnetic mirrors. Third, in the strongly coUisional limit: A <C Pe, the particles do not move along the 
field lines, magnetic mirroring is inefficient and, thus, kb approaches the classical Spitzer value. Note also that in the first 
and second cases, the heat flux is anisotropic for small-scale perturbations, kls 2> 1, because the magnetic field is locally 
ordered. The flux is given by equation (ph with k = kb across the fleld lines, and by equation (O) or equation (pf) with 
1^ = i^Sp, depending on plasma coUisionality, along the field lines. 

Now we estimate the fraction of particles, 'd, which can pass through magnetic mirrors. Let us consider a particle which 
moves along a field line from a region "1" with a weak field to a region "2" with a stronger field. If the field gradient scale 
is much larger than the Larmor radius (which is the case for most astrophysical systems), the particle's magnetic moment 
/i = mv\/2B {y± is the velocity component perpendicular to the local field) is conserved as an adiabatic invariant. That 
is mv\^/2Bi = 771^^2/2-82- The energy of the particle is also constant. At the point where the particle is refiected its 
parallel velocity vanishes, therefore v\2 — '^\i +^yi- Combining these two equations we obtain that particles are refiected 
from the magnetic mirror if their pitch angles are greater than the minimal pitch angle, defined as: 

Sin^ 6I,„ = Smin/Smax- (A4) 

The particles with 6 < 6m, which can pass though the mirror, form a "loss cone", the solid angle of which is 

flic = 27r(l — COS^ni). 

There is a second loss cone in the opposite direction along the field line, due to symmetry. Assuming an isotropic 
distribution of particle velocities (in the weak field region), we estimate the fraction of free streaming particles as 

d = 2nijAn = 1 - (1 - B„,i„/S„,ax)'/' ~ 0.3, (A5) 

where we assumed that the amplitude of magnetic turbulent fiuctuations in an accretion flow is typically SB ^ (B), i.e., 



In hot accretion flows, the Coulomb mean-free-path is often much larger than the size of the flow, so that A > Lrji. It 
is convenient then to write the conductive heat flux as in equation (|l4|). 



ff, w>»e „..^j:„,2;g^....^,.,,„ ,M, 



^^2 ^.2 
'J'cond = —a. 

iiK ax 

i—e^p 



is the coefficient of thermal conduction, i — e,p because both species contribute to heat transport, and F{e,p) = 
m.ax [F (e) , F (p)] properly corrects for the two-temperature flow, and where H/R ~ Cs/vg — Cs/{V^^kR) is of order 
unity in hot flows, H is the accretion disk scale height and wg is the free-fall speed. We also used here that vxe.p — Cse,p- 
Thermal conduction is dominated by electrons in a one-temperature flow because they are lighter than protons. In a two- 
temperature flow, however, the protons may be much hotter than the electrons and, hence, will dominate the conduction. 
We write the relative contributions of electrons and protons to thermal transport in two- and one-temperature zones: 

rrie Cse dc^^/dx ^ ,,,„e 

nip Cs dcl/dx ' ^ 2\/2 V mp ~ ' ' z^^-, 

F{,p)2T = 1, 




where c^ = c^^ -f {me/mp)c'^^ = nkB{Tp + T^). We have made use of the fact that in a two-temperature zone (Tp 3> T,.) 
Cj ~ c^ and in a one-temperature zone {Tp — Te) c^ — 2c^ . 

PROPERTIES OF THE HOT SETTLING FLOW ONTO A NEUTRON STAR 



The flrst ana lytical solution for a hot accretion flow settling onto a rapidly rotating NS was found by Medvedev fc 



Narayan (2001 ). We briefly present here some relevant results and derive some new results. The hot settling flow is 
geometrically thick with the vertical thickness being comparable to the local radius of the flow, H ^ R. The flow is hot 
with the protons being at nearly the virial temperature. The flow rotates with a sub-Keplerian angular velocity. The 
flow is optically thin and is cooled by free-free emission. The flow is powered by the rotational energy of the central star; 
hence the flow parameters (all but the radial velocity) are independent of the mass accretion rate. As in the case of an 
ADAF, the mass accretion rate must not exceed a critical value which is of order a few percent of the Eddington value; 
otherwise, the cooling is too strong and a thin Shakura-Sunyaev disk forms. 

The settling flow consists of two zones: an inner two-temperature zone for r <^ 10^'^ (r — R/Rs, where Rs = 
2.95 X lO'^m cm is the Schwarzchild radius), and an outer one-temperature zone for r ^ 10^'^. Each zone is described by 
a self-similar solution. The two- and one-temperature solutions read: 

P2T — 8.9 X 10~^TO~^as^r~^g/cm^, piT ~ 4.9 x 10~^TO~^as^r~^ g/cm^, 

Op,2T = (l/6)r-l, 0p,lT = 9p,2T/2, 

^e,2T ^ llr~l/2, ee,lT ^ ISSr"!, (Bl) 

Q2T ^ 7.2 X lO^m-^sr'^/^ rad/s, Qit = ^2T, 

V2T — 1-4 X \{f'rna^^s^'^r^ cm/s, vit — 2.6 x \{fma~^s~^r'^ cm/s, 

where Op = ksTp/mpC^ and similarly for the electrons, m — Mns/Mq, rh — M / M-^dd is the mass accretion rate in 

Eddington units, MEdd = 1-4 x lO^^m g/s, a is the viscosity parameter, and s = Qns/^k{Rns) is the NS spin in units of 
the Keplerian angular velocity at the NS radius. Note that, except for v, none of the quantities depend on the accretion 
rate to. The one-temperature self-similar solution is valid from r ^ 10^-^ to approximately 

r,s -- {TAa'^s'^/mY - 5.5 x lO^ai^s^rhzl- (B2) 

Beyond this radius the solution is not self-similar. 

The above self-similar solutions were obtained without taking into account thermal conductivity of the gas. As we 
will now see, thermal conduction in a hot settling accretion flow is very large, so that the thermal flux is energetically 
important and modifles the structure of the cqulibrium settling flow itself. 

To be consistent with our simple analysis of the thermal instability in §y, it is sufficient to consider the one-temperature 
regime only. The hot settling solution is obtained from the condition that the total local heating rate is equal to the total 
cooling rate. Including the contribution from the divergence of the conductive flux, the condition now reads 

where Qcond, Q-, and Q^ are given by equations (|l4|), (^. Since the additional hea t flu x term has exactly the same 



radial dependence as the other two terms, the power-law scalings given by equations (Bl) remain valid, but with their 



pre- factors change d. Let us write the new solution as p 
balance equation ([B3|) takes the form: 



pr 



Op ~ OpT ^, etc.. With these definitions, the energy 



:Pc 



l^KT 



'9/2 



as 



+ UpaA^Cp\el/\-''\ 



^' ' '^ V 3 

Clearly the effect of thermal conduction is to to re-define the quantity as^ to 



(B4) 



as 



as 
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(B5) 



where we have used Op = 1/6. With this re-definition, we may continue to use equations (Bl) for the equilibrium flow. 

Next, we estimate the typical value of the thermal conduction coefficient in the hot settling flow. From equations (15) 
and ( |A7| ) we have 

R 



2-1/2 _i^^^(e,p) 



10-2^_it?_i^(e,p). (B6) 

From (A7) and (Bl) we obtain F{e)2T — O.lSr'^/^ < 10.9 because the two-temperature zone exists for r < 10^'^. That is, 
in the two-temperature zone F{e,p) changes from ^^ 1 in the inner parts to ~ 10 in the outer parts and then smoothly 
approaches the value of ~ 15 in the one-temperature zone. Thus, the overall quantity F{e,p) = max.[F{e), F{p)] lies 
in the range 1 < F{e,p) <; 15 for the entire flow. Moreover, we note that F{e) <, F{p) only for small radii r ^ 13, 
where the self-similar settling solution is not terribly accurate because of Comptonization. Therefore, we see that thermal 
conduction is dominated by the electrons over nearly the entire accreti on fl ow. 

Finally, we demonstrate that (|l^) is the relevant limit of equation ( |A3| ) and, hence, that the above estimate for the 
thermoconduction coefficient is accurate. For this, we must demonstrate that A ^ L^jji in the settling accretion flow (see 
equation A3). 

From the estimates above we have as^ + {2/9)ac ~ 2.2 x 10^^S,-i'&~iF{e,p). From equation (pi]), the electron mean- 
free-path, A ^ 0.579^/ p, normalized by the local radius of the flow, is estimated in both zones to 



X2t/R ^ 12 [e-ii9-iF(e,p)]-' , Xit/R ^ 17 [e-i^-iF(e,p) r^]'^ 



(B7) 



The proton mean free path is comparable to or larger than the above estimate. To calculate Lm^ the Larmor radius of 
the electrons is needed. It may be estimated as follows. The magnetic pressure Pm = Pgas/P yields 



B ~ {8npcl/P) 



1/2 



Since c^ — c^^ + {melrap)c^g^ ^ c^p to within a factor of two and the thermal velosity of particles is VTe 



PL ^ VTe/^B ^ 1.1 X 10-« ( ^ 



1/2 



we obtain 



(B8) 



wher e 17 b — e B/rripC is the electron cycloron frequency. For a typical /3 ~ 10 (see e.g., Medvedev 2000 ; Quataert fc 
Gruzi lov 1999| ) , the normalized electron Larmor radii in the two- and one-temperature zones are: 



Pl,2t/R - 2.2 X 10-iV-i/2ri/4 [^_,^_,F{e,p)] 



-1/2 



Pl,it/R ' 



5.1 X 10-i"m-i/2 [^_^§_^Fie,p)] 



-1/2 



(B9) 



Since the two quantities differ at most by a factor of order unity, we use Pl/R ~ 5 x 10 i" as a representative value. Now 
we can estimate the Rechester-Rosenbluth length from equation (|A2|), 



Finally, we have 



Lnn/R^^lniC/iPL/R)] - 10-1^-1(19 + In^-i) ^ 1.9^_i. 



X2T/LRR~6.3[Cid^iF{e,p)] \ \it/Lrr ^ 8.9 [C\^-iF{e,p)] ' 



(BIO) 
(Bll) 



We see that in the two-temperature zone, the mean free path is large, A ^ Lrr, so that equations (g_3|) and ( |A3| ) are 
justified. In the one-temperature zone, the result is marginal: A ~ Lrr. However, considering that (i) the transition 
between the regimes with A <C Lrr and A 3> Lrr is not sharp in reality, (ii) the "boundary" itself, A ~ Lrr , is uncertain 
to within a numerical factor of a few, and (iii) the one-temperature zone is rather small [c.f., equation (B2)], we believe 
that we may use equation ( p^ ) throughout the flow. The error will not be greater than a factor of order unity. 
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PROPERTIES OF THE SLE SOLUTION 

The SLE slim disk solution was the first hot accretion solution found. The SLE solution is optically thin. It is cooling- 
dominated via free- free emission (the original paper by SLE considered Compton cooling). As in other hot solutions, 
there are both two- and one-temperature zones in this flow. In the two-temperature zone, the normalized electron and 
proton temperatures are related as follows: 



TrlnA 



12a 



1/2 



e]!'^ ~ 230y^ 



(CI) 



where a/ is the fine structure constant and In A — 15. This equation is easily obtained from the equality of the electron 
cooling rate via Bremsstrahlung and the energy transfer r ate from the protons to the e lectrons via Coulomb collisions. 



T„ 



Qcoui; the expressions for these r ates are given in Medvedev fc Narayan (200l| ). In the one-temperature zone, 
{mp/me)Op. Following Shapiro, Lightman, fc Eardley (1976|), we obtain the following solutions: 



= Te, i.e 

Cs,2T 

P2T 

iH/R)2T 



6.5 X lO^m^/^'a ^^^r ^/^ cm/s, 

2.9 X lO-^TO-imi/^ai/^r-is/s g/cm^ 

3.1 X 10-imi/4a-i/2j,i/8 



Cs,lT 

Pit 
{H/R)iT 



1.9 X 109mi/4a-i/2^-3/8 (,jj^/g^ 



1.2 X IQ-'^m-^m^/^a^/'^i 
2.2 X 10-imi/4a-i/2^i/8 



-15/8 



/ 3 

g/cm , 



(^2) 
From the above equations, the dimensionless temperature in the two-temperature zone (r <J 10^) is 9p = [cs/cf' ~ 

0.05 r~^/^ which is significantly smaller than the temperature of the hot settling flow in which dp ~ 0.2r^^. In the 
one-temperature zone, thetap ~ 0.004 r~^/^ is again smaller that the temperature of the hot settling flow, 6p ~ 0.1 r^^, 
up to a very large radius r ~ 4 x 10^. We see that the temperature of the gas in the SLE disk is smaller than in the hot 
settling flow. Therefore, we expect that collisions are more frequent and they may change the thermal properties of the 
gas. Indeed, we denonstrate now that A <^ Lfm in the SLE disk, that is, thermal conduction is in the semi-coUisionless 
regime (see Appendix ^ . 

From equations (O) with T = c^irip/kB and (CI), we obtain the normalized coUisional mean-free-paths 



^2T.p/R — 2.9m_2 a_i 



3/4^-5/2^-5/8 



^2T,e/R ■ 



0.3m_2 Q^-1 ^^2 



AiT/i?=i2.2x 10 



-2 • 3/4 -5/2 -5/8 



The normalized electron Larmor radius follows from equation (B8): 

PL.2T/R ^ 1-3 X 10-'^°m-^/^mZl^^rl^^, Pl.it/R ^ 3.6 x lO^^m" 



1/2 • -1/8 -1/4 -1/16 



(C3) 



(C4) 



where we take r ~ 100 as a characteristic radius at which a one-temperature flow becomes two-temperature. Since 
pL oc (Tto)^/^, the proton Larmor radius is ^mpjmf, ^ 43 times larger. We estimate L_r/j from equation ( |A2|) as 



iW^~-ein(pL/i?)~2.5C_i 



(C5) 



for the electrons and it is slightly less for the protons. Comparing this with (C3), we see that A ^ i_R_R in that part of the 
disk where the electron thermo-conduction dominates. Only in the inner parts of the disk, where the proton conductivity 
becomes large, do we have A J> L^fl. 



Finally, we calculate the thermal conductivity of the gas in the SLE disk. It is determined by equation (A3) with 
A < ifiij: 

HB = Ksp^ils/LRR) = Ksp^/I Hpl/R)\ ~ 4 X IQ-^^^iKsp < Ksp- (C6) 

This thermoconductivity is about two orders of magnitude smaller than the Spitzer value. As we discussed above, the 
protons may significantly contribute to the thermal conduction flux and may even dominate the electron contribution 
in the inner parts of the disk. The proton contribution, which is nspp = y/me/mpKspeiTp/Te)^^^ (see equations pffl, 
pT[|), becomes dominant in the two-temperature zone with Tp ^ Tg. In o ne-t emperature zone, the electron contribution 
dominates, since Kspp — (l/43)Kspe- Using the self-similar SLE solution (C2), we calculate the nominal Spitzer thermal 
conductivity in the SLE disk 



/.Sp,2T ^ 2.7 X 102im5/4a-5/2r-i5/8^ 



«;Sp,iT ^ 2.5 X 102"m5/4a-5/2^-i5/8 



(C7) 



in CGS units. The real value of k is less than the Spitzer value by the factor given in (C6). 
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Fig. 1. — Real part (left panel) and imaginary part (right panel) of the frequencies of three modes. The curves labeled 1 and 2 refer to the 
two acoustic modes, and the curves labeled 3 refer to the thermal mode. The dashed curves correspond to the dispersion relation (El), which 



does not include thermal conduction. The following parameter values were used: Kepi/f^if 



/0.5, 



, Tcool.ff 



ion: 



The mode frequencies are normalized by the Keplerian frequency. The solid curves correspond to the dispersion relation 
thermal conduction. Here, r^ond = """cooLffi and the other parameters are the same as before. 



Ca = UkR, 7 = 5/3. 
iVw, which includes 





